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I. INTRODUCTION 

The interest in luminescent materials has increased 
in recent years because of their numerous technological 
applications.El Most of these materials involve the doping 
of a pure ionic crystal, that is substitution of some of 
the ions by other ions with specific absorption-emission 
characteristics (Tl + , Ce 3+ , etc). The presence of impu- 
rity ions in ionic crystals induces geometrical distortions 
in the host lattice. Those distortions are an important 
ingredient in the scintillation process and thus it is im- 
portant to establish a reliable description of them. 

A convenient way to calculate the lattice distortions is 
to model the doped crystal by a finite cluster centered 
on the impurity and embedded in a field representing 
the rest of the host lattice. This cluster approach has 
been used to studv^the geometrical and optical properties 
of doped crystals.ErEj The cluster (active space) can be 
studied by using standard quantum-mechanical methods. 
The rest of the crystal (environment) can be described in 
several ways. In the simplest and most frequently used 
approach, the environment is simulated by placing point 
charges on the lattice sites, but this procedure has to be 
improved in order to obtain a realistic deapaptiqa of t,t 
lattice distortions around the impurity.tTBEJoOl 
Model potentials have been developed to represent the 
effects of the environment on the active cluster, that in- 
clude attractive and repulsive quantum-mechanical terms 
aside from the classical Madelung term,t3 but a problem 
still remains: the large computational cost of conven- 
tional molecular orbital (MO) calculations prevents from 
performing an exhaustive geometrical relaxation of the 
lattice around the impurity. In the most accurate MO 



calculations, only the positions of the ions in 

the first shell around the impurity are allowed to re- 
lax. However, geometrical relaxations far beyond the 
first shell of neighbors can be-exeected. In fact, recent 
classical simulations of solidsJ.TP.i performed employing 
phenomenological potentialSjEHcll have shown the impor- 
tance of considering appropriate large-scale lattice relax- 
ations in the study of a variety of intrinsic and extrinsic 
defects in ionic crystals. 

In this contribution we report theoretical calculations 
of the lattice distortions induced by a Tl + impurity in 
two well-known ionic scintillators, namely NaI:Tl + and 
CsI:Tl + . For this purpose we use the ab initio Per- 
turbed Ion (PI) model,E3~Ej which circunvents the prob- 
lems mentioned above: (a) The active cluster is em- 
bedded in an environment represented by the ab initio 
model potentials of Huzinaga et al.cB; (b) The computa- 
tional simplicity of the PI model allows for the geomet- 
rical relaxation of several coordination shells around the 
impurity.tifli-3 

The remainder of this paper is organized as follows: In 
Section [| we describe the different clusters which have 
been used in order to assert the influence of the embed- 
ding scheme and of the size of the active cluster on the 
results. In Section 111 we present and discuss the results 
of the calculations, and Section IV summarizes the main 
conclusions. 
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II. CALCULATIONAL STRATEGY. 
DESCRIPTION OF THE DIFFERENT ACTIVE 
CLUSTERS 

The ab initio Perturbed Ion model is a particular appli- 
cation of the theora of electronic separability of Huzinaga 
and coworkers J£3c3 in which the basic building block is 
reduced to a single ipa. The PI model was first developed 
for perfect crystals Its application to the study of im- 
purity centers in ionic crystals has been described in refs. 
[ll],|l4], and we refer to those papers for a full account of 
the method. In this case, an active cluster containing the 
impurity is considered, and the Hartree-Fock-Roothaan 
(HFR) equationsCJ of each ion in the active cluster are 
solved in the field of the other ions. The Fock operator in- 
cludes, apart from the usual intra-atomic terms, an accu- 
rate, quantum-mechanical crystal potential and a lattice 
projection operator which accounts for the energy con- 
tribution due to the overlap between the wave functions 
of the ionsrJ The atomic-like HFR solutions are used to 
describe the active cluster ions in an iterative stepwise 
procedure. Also, as we shall see later in more detail, the 
positions of some ions in the active cluster are allowed 
to relax. The wave functions of the lattice ions external 
to the active cluster are taken from a PI calculation for 
the perfect crystal and kept fixed during the embedded- 
cluster calculation. Those wave functions are explicitely 
considered up to a distance d from the center of the ac- 
tive cluster such that the quantum contribution from the 
most distant frozen shell to the effective cluster energy is 
less than 10~ 6 hartree. In practice, this means 25 shells 
for NaI:Tl+ and 21 shells for CsI:Tl+. Ions beyond d con- 
tribute to the effective energy of the active cluster just 
through the long-range Madelung interaction, so they are 
represented by point charges. At the end of the calcula- 
tion, the ionic wave functions are selfconsistent within the 
active cluster and consistent with the frozen description 
of the rest of the lattice. The intraatomic Coulomb corre- 
lation correction, which is neglected at the Hartree-Fock 
level, is computed by using-nthc Coulomb-Hartree-Fock 
(CHF) model of Clementi.ElE3 An effective energy can 
be assigned to each ion that collects all the contributions 
to the total crystal energy involving that ion.ca 

In this work we have employed several active clusters 
of increasing complexity and with different embedding 
schemes. Now we describe each cluster (see fig. 1). First 
for NaI:Tl + , that is, a substitutional Tl + impurity in a 
Nal crystal with the rock-salt structure: 

• A) (T1I 6 ) 5 ~:PC. This active cluster is formed by 
the impurity cation Tl + and its 6 nearest neigh- 
bors, whose positions are allowed to relax. The rest 
of the crystal is simulated by using point charges 
(PC). Notice that in this case, and also in clus- 
ter A* below, the environment is simplified with 
respect to the statements at the beginning of this 
section, in the sense that the "frozen ions" are just 
point charges. 



• B) (T1I 6 ) 5 ":NaI. The active cluster is the same as 
in A, but now the environment is represented by 
ions whose "frozen" wave functions are obtained 
from a PI calculation performed for the pure Nal 
crystal. In the perfect crystal all the cations (or 
anions) are equivalent, and there is a single nearest- 
neighbor Na-I distance. These characteristics are 
lost around the substitutional impurity. 

• C) (TlIi 4 Nai 8 ) 5+ :NaI. The active cluster is formed 
by 33 ions corresponding to the central impurity 
cation plus four coordination shells. The rest of the 
crystal is represented by frozen ions with PI wave 
functions as in B. In the calculations performed on 
this cluster only the positions of ions in the first 
shell around Tl + are allowed to relax. The other 
three shells in the cluster provide an interface be- 
tween the first shell and the frozen environment. 
The ions in this interface respond to the distor- 
tion induced by the impurity by adapting selfcon- 
sistently their wave functions to the new potential, 
but not their positions. 

• D) (TlI 92 Na86) 5 ~:NaI. This is the largest active 
cluster used for this material. It has 179 ions which 
correspond to the central Tl + cation plus 12 coor- 
dination shells. Only the 33 ions of the four inner 
shells relax their positions; the eight outer shells 
provide the interface between the inner ones and 
the environment of frozen PI ions. 

Next, for CsI:Tl + , that is, a substitutional Tl + impurity 
in a Csl crystal with the CsCl structure, we define four 
clusters in a similar way: 

• A*) (T1I 8 ) 7 ~:PC. The active cluster is formed by 
the impurity plus its eight nearest neighbors, and 
is embedded in an environment of point charges. 

• B*) (Tils) 7 ~:CsI. This cluster is similar to A*, but 
the environment is now simulated by frozen ions 
with wave functions taken from a PI calculation 
for the pure Csl crystal. 

• C*) (Tll3 2 Cs 32 ) +:CsI. The active cluster has 65 
ions and is formed by the impurity plus six coor- 
dination shells. Only the positions of the ions in 
the first shell around Tl + are allowed to relax. The 
environment is the same as in B*. 

• D*) (TlI 88 Cs 92 ) 4+ :CsI. The active cluster has 181 
ions: the impurity plus twelve coordination shells. 
The ions in the six innermost shells can relax their 
positions, and the environment is as in B*. 

The geometrical relaxations of the shells in the differ- 
ent cluster models described above have been performed 
by allowing for independent breathing displacements for 
each shell of ions, and minimizing the total energy with 
respect to those displacements until the effective cluster 
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energies ace converged up to 1 meV. A downhill simplex 
algorithmEj was used. For the iqes we have used large 
STO basis sets: (5s4p) for Na+,Eil (lls9p5dUjor 
(Ils9p5d) for Cs+,El and (13s9p7d3f) for T1+.E2 

III. RESULTS AND DISCUSSION 

Before presenting the results for the lattice distortions 
around the Tl + impurity we test the consistency of the 
embedding method for the case of pure (undoped) crys- 
tals. For this purpose we compare the results of two 
sets of calculations for Nal. In the first one we study the 
clusters labelled A, B, C, D in Section [n[ with the em- 
bedding scheme indicated there for each case, but with 
a Na + cation instead of the Tl + impurity. This is, of 
course, just the case of the pure Nal crystal treated by 
the embedded cluster method. One can then compare 
the results with those of a standard PI calculation for 
the perfect crystal.c3 We can anticipate some differences 
between the two methods since in the usual PI calcu- 
lation for a perfect crystal all the cations (or anions) 
are equivalent, while in the embedded-cluster descrip- 
tion of the same system the cation acting as a fictitious 
impurity and the other cations of the crystal are not de- 
scribed in the same way. This systematic error is, in fact, 
what we want to remove in the analysis of the distor- 
tions around the true, Tl + , impurities. If we call Ri the 
distance between the Na + cation and its first I - neigh- 
bors, Table I gives Ri for the four cluster models A, B, 
C, D, together with the difference ARi = [Ri (cluster) 
- Rf 1 (crystal)] between the embedded-cluster result and 
that of the perfect crystal at equilibrium and the relative 
deviation Ai?i/Rf 1 (crystal). The PI model predicts for 
the crystal an equilibrium value Rf 7 (crystal) of 3.237 

A, in very good agreement with the experimental value 
pjexp _ g 240 A. The first two cluster schemes, A and 

B, give large distortions AR%. In contrast, a contrac- 
tion of Ri smaller than 1 % occurs for cluster models C 
and D. The conclusions for Csl are similar. A good self- 
embedding is only achieved for clusters C and D (C* and 
D*). Comparison of B and C shows the importance of 
a smooth interface between the inner cluster core, where 
ions are allowed to move, and the frozen environment 
around the cluster. Comparison of C and D establishes 
the necessity of allowing for relaxation of the radii of 
several shells around the impurity. The systematic error 
that the embedding-cluster scheme makes has to be taken 
into account in order to interpret properly the distortions 
induced by the Tl + impurity. 

Now we present the results for the lattice distortions 
induced by the Tl + impurity in Nal and Csl. In this case 
we have calculated 

AR, = R,{NaI : Tl + ) - R z (NaI : Na + ), (1) 

where Ri (i=l, 2, 3, ...) refer to the radii of the first, sec- 
ond, ... shells around the impurity and both calculations 



have been performed in the embedded cluster scheme. 
In this way the systematic errors of the cluster method 
analyzed in Table I tend to cancel. As before, Ri is the 
distance between the central cluster cation, Tl + or Na + , 
and the I - anions in its first neighbor shell. For the rea- 
sons given above we only trust the results from clusters 
C, D (or C*, D*) and the results for Ri and ARi are 
given in Table II, which also contains results for Tl + in 
Csl. Although we discourage the use of models A and 
B, nWst for the purposes of comparison with Berrondo et 
al.t3 we have also calculated Ai?i for NaI:Tl + in model 
A and for CsI:Tl + in model A*. Berrondo used a simi- 
lar embedded-cluster scheme, although his calculational 
method was different, based on the linear combination of 
atomic orbitals (LCAO). Taking as reference the exper- 
imental nearest-neighbor Na + -I~ distance in the bulk, 
Berrondo obtained ARi — 0.54 A, and taking the same 
reference we obtain AR 1 = 0.50 A. The good agreement 
between the two calculations provides a check of our theo- 
retical method. However this large value of ARi is partly 
due to the fact that the calculated Ri(NaI:Tl + ) and the 
reference R x (pure crystal) are not fully consistent with 
each other for the reasons discussed above. If we cal- 
culate Ri(pure Nal) also by using model A (with Na + 
instead of Tl + ), that is, like in eqn. |l|, then a corrected 
"smaller" value AR 1 = 0.20 A (or 5.6 %) results. 

Returning to models C and D we predict, using eq. ^ 
a lower expansion of Ri. In fact, Ai?i is equal to 0.039 
A (1.2 % expansion) for model D. However, the distor- 
tions are not restricted to the first shell. The second and 
fourth shells also expand a little, while the third shell 
suffers a contraction. The expansion of the first and sec- 
ond shells is consistent with the fact that the ionic radius 
of Tl + , equal to 1.40 A, is larger than the ionic radius 
of Na + (0.95 A). In summary we find that the struc- 
tural relaxation around the Tl + impurity in Nal is small, 
although more complicated that usually assumed by pre- 
vious workers. In fact our calculations suggest that the 
relaxation could even go beyond the fourth shell. 

Let us now turn to CsI:Tl + . The calculation with 
model A* would give a small expansion ARi of the first 
shell if we take as reference R^ p (pure crystal), and our 
calculation would, then be in agreement with that of 
Berrondo et al.!E3 However the ionic radius of Tl + is 
lower than that for Cs + , so an expansion of Ri is not 
to be expected. By using, instead, eq. |l| we obtained a 
corrected result ARi = —0.062 A, that is a contraction. 
Cluster model C*, in which only the first shell around 
Tl + is allowed to distort, gives a small expansion of Ri 
and finally model D* recovers again the expected con- 
traction. Beyond the first shell, ARi oscillates (Ai?2 and 
Ai?4 are also negative, while AR 3 , AR 5 and AR 6 are 
positive), although the values of ARi are very small, in 
fact smaller than in NaI:Tl + . 

The conclusion from our calculations is that the distor- 
tions around the impurity affect ions in several coordina- 
tion shells around the impurity, and that consideration 
of both the atomic and the electronic relaxation of ions 



3 



in those shells is required. 

For each active cluster we have also calculated the ab- 
sorption energy corresponding to the intra-atomic tran- 
sition from the singlet ground state to the triplet excited 
state of the thallium ion (6s 2 ( 1 S) -» 6s 1 6p 1 ( 3 P)). Fol- 
lowing the Franck-Condon principle, this absorption has 
to be calculated at the frozen ground state geometrical 
configuration. Thus, to obtain the energy of— the local- 
ized excitation, we solve the HFR equationsca for that 
geometrical configuration of the active cluster, with the 
T1+ ion in the excited electronic state 6s 1 6p 1 ( 3 P), and 
evaluate the absorption energy as a difference of effective 
ionic energies £3 

E abs = E eff [Tl+( 3 P)] - E eff [Tl+CS)}. (2) 

To calculate the contribution of an electronic open shell 
to the intraatomic energy of a LS electronic configura- 
tion, the Hartree-Fock-Roothaan formalism requires the 
coupling constants of that specific LS term as input. For 
the description of the open shell ion we have used the 
coupling constants of the 3 P term as given in ref. |5^ and 
the same basis set as for Tl+(6s2). The absorption spec- 
tra of these doped crystals are more structured (there 
are several absorption bands); our calculated transitions 
should be identified with the band A of the absorption 
spectrum of NaI:Tl + (4.25 eV),E3 and. with the first ab- 
sorption band of CsI:Tl+ (4.27 eV)H Both in NaI:Tl+ 
and CsI:Tl + , the absorption energies we obtain become 
closer to the experiment as we improve the description 
of the local geometry around the impurity, as Table II 
shows. In particular the absorption energies obtained 
with models D and D* are remarkably accurate, and this 
can be ascribed to the accurate representation achieved 
for the lattice distortion around Tl + . Treatment of other 
absorption bands would require to depart from some of 
the basic assumptions of the PI model. 



the central impurity plus its first shell of neighboring ions 
only, should be taken with some care. In CsI:Tl + , our 
results indicate that only the first shell suffers an appre- 
ciable displacement. Calculated absorption energies are 
in excellent agreement with the experimental values, and 
that agreement improves as the description of the cluster 
becomes more realistic. 

Thus, the PI model is able to give a realistic descrip- 
tion of the geometrical distortion around a defect in an 
alkali halide crystal. Of course, the advantages come with 
some drawbacks: The simplifying assumptions present 
in the actual version of the PI model make it nonsuit- 
able for the study of excited states involving interatomic 
charge transfer from the impurity. The method also has 
difficulties for accounting for the crystal-field splitting 
of open shell excited states. These drawbacks concern- 
ing the treatment of excited states should not overlook 
the essential message of this paper, that is, that the PI 
method combined with an embedded cluster model allows 
to obtain accurately the distortion of the lattice around 
impurity centers in ionic crystals. A combined strategy, 
namely, PI calculations to obtain the distortions, followed 
by a MO calculation on the optimized geometry, can be 
a desirable option in the future. 
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IV. SUMMARY 



We have reported a study of the local lattice distortions 
induced by a Tl + impurity in NaI:Tl + and CsI:Tl + scin- 
tillators. To that end, the ab initio Perturbed Ion (PI) 
model has been used. Large active clusters, embedded in 
accurate quantum environments representing the rest of 
the crystal, have been studied. The importance of per- 
forming parallel calculations for the pure systems in order 
to suppress systematic errors from the calculated local 
distortions has been stressed. The self-embedding per- 
formance of the method has been also analyzed with the 
result that it is necessary to include in the active cluster 
some shells of ions with frozen positions but able to adapt 
selfconsistently their wave functions to the displacements 
of the inner shells. The local distortions obtained extend 
beyond the first shell of neighbors in NaI:Tl + . Thus, the 
assumptions frequently employed in impurity-cluster cal- 
culations, which consider the active space as formed by 
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Captions of Tables 

Table I. Self-embedding test. Values of Ri, the ra- 
dius of the first coordination shell around Na+ in Nal 
obtained from embedded-cluster calculations, and their 
deviations ARi from the theoretical value obtained in a 
PI calculation for the perfect crystal. Relative deviations 
ARi/Rf 1 (crystal) are given in brackets. Note the good 
selfembedding achieved in cluster models C and D. 

Table II. Distortion induced by a Tl + impurity in 
Nal and Csl. Distances R^ (i=l, 2, ...) in A between 
the Tl + and its succesive coordination shells are given, 
as well as the differences AR^ with respect to the values 
obtained in calculations for the perfect crystal (also with 
the cluster model) and percent distortions. Last column 
gives absorption energies (see text) together with percent 
deviations from the experimental value. 



Captions of Figures 

Figure 1. The different active clusters studied. 
Small spheres represent cations and large spheres rep- 
resent anions, (a) (T1I 6 ) 5 ~: (b) (TlIi 4 Nai 8 ) 5+ ; (c) 
(TlI 92 Na 86 ) 5 "; (d) (T1I 8 ) 7 ~; (c) (TlI 32 Cs 32 ) +; (f) 
(TII88CS92) 4+ . These clusters are embedded in the crys- 
tal as indicated in the text. 
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NaI:Na+ Cluster Ri (A) ARi (A) 



A 


3.540 


+ 0.303 






(+ 9.36 %) 


B 


3.159 


- 0.078 






(- 2.41 %) 


C 


3.213 


- 0.024 






(- 0.74 %) 


D 


3.234 


- 0.003 






(- 0.09 %) 



TABLE I. 



G 



[ 



NaI:Tl+ 


Cluster 


Ri R 2 


R 3 


R 4 ARi 


AR 2 


AR 3 


AR 4 


Abs. (eV) 




A 
C 
D 


3.738 
3.274 

3.273 4.692 


5.510 


+ 0.198 
(+5.59 %) 

+ 0.061 
(+1.87 %) 
6.612 + 0.039 
(+1.21 %) 


+ 0.075 
(+1.62 %) 


- 0.054 
(-0.97 %) 


+ 0.043 
(+0.64 %) 


3.91 (-8.00 %) 
4.19 (-1.41 %) 
4.23 (-0.47 %) 




CsI:Tl+ 


Cluster 


ARi 


AR 2 


AR 3 


AR 4 


AR 5 


AR 6 


Abs. (eV) 




A* 
C* 
D* 


- 0.062 
(- 1.62 %) 

+ 0.013 
(+ 0.30 %) 

- 0.047 
(-1.14 %) 


- 0.008 
(- 0.15 %) 


+ 0.003 
(+ 0.06 %) (- 


- 0.016 + 0.006 
0.21 %) (+ 0.09 %) 


+ 0.001 
(+ 0.01 %) 


4.17 (-2.34 %) 
4.21 (-1.40 %) 
4.26 (- 0.23 %) 



TABLE II. 
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